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Multicritical behaviour of the compressible systems 
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The behaviour of uniform elastically isotropic compressible systems in critical and tricritical points 
is described in field-theoretical terms. Renormalizationgroup equations are analyzed for the case 
of three-dimensional systems in a two-loop approximation. Fixed points corresponding to various 
types of critical and multicritical behaviour under various macroscopic conditions imposed on the 
system are distinguished. It is shown that the effect of the elastic deformations on the critical 
behaviour of compressible systems is significant. It manifests itself both in a change in the critical 
exponents of Ising magnetic and in the appearance of multicritical points in phase diagrams at any 
dimension of the order parameter. It is also shown that, in a number of experimental investigations, 
the multicritical behaviour is not tricritical, as it has been stated, but tetracritical. The influence 
exerted by elastic deformations on systems with phase diagrams already containing multicritical 
points is analyzed. It is shown that the effect of elastic deformations leads to a change from bicritical 
behaviour to a tetracritical one. 

PACS numbers: 64.60.-i 



I. INTRODUCTION 



Phase diagrams of many materials exhibit multicrit- 
ical behaviour when several lines of phase transitions 
are crossed. The tricritical points, in which the second- 
order phase transition is replaced by the first-order tran- 
sition, are of particular interest. Experimental studies 
are mostly concerned with systems subjected to exter- 
nal pressure, i.e. those influenced by elastic deforma- 
tions. Owing to this circumstance, taking into account 
the relationship between the order parameter and elas- 
tic deformations is important. As first shown in [1], in 
the elastically isotropic case, the critical behaviour of 
compressible systems with quadratic striction is unsta- 
ble with respect to the relationship between the order 
parameter and acoustic modes, and a first-order phase 
transition close to the second-order transition occurs in 
systems of this kind. However, the conclusions of [1] are 
only valid at low pressures. Itwas shown in [2] that at 
high pressures, a more fundamental influence is exerted 
on the system by deformational effects induced by the ex- 
ternal pressure, beginning from a certain tricritical value 
P t . This changes the sign of the effective constant of 
interaction between fluctuations of the order parameter, 
and, as a consequence, leads to a change of the type of 
the phase transition. For uniform compressible systems, 
two types of tricritical behaviour and the existence of a 
tetracritical point, at which two tricritical curves meet, 
were predicted in [2]. Calculations made in terms of a 
two-loop approximation [3] confirm the existence of two 
types of tricritical behaviour for the Ising systems and 
yield the tricritical exponents. 

According to the criterion obtained in [1], striction 
effects, regarded as additional thermodynamic parame- 



ters, change the mode of the critical behaviour only in 
systems with a singular behaviour of heat capacity in 
the absence of deformations. The specific heat exponent 
a(C ~ \T — T c \~/ alpha ) is only positive for Ising-like mag- 
netic. In the XY-model and the Heisenberg model, the 
specific heat exponent is positive for "rigid" systems, and, 
therefore, elastic deformations must have no effect on the 
critical behaviour. Hence follows that the critical value 
of the order parameter dimensionality n c < 2. 

In this study, we develop further the model of phase 
transformations in uniform compressible systems char- 
acterized by various dimensions of the fluctuating order 
parameter [4, 5] . We analyze these transformations using 
the renormalization-group methods in a two-loop approx- 
imation, directly in the three-dimensional space. Also, 
we consider the conditions for occurrence of a tricritical 
behaviour owing to effects of long-range interaction of the 
order-parameter fluctuations caused by long-wavelength 
acoustic modes. Since the dependence of the exchange 
integral on distance makes a major contribution to the 
striction effects in the critical region, we consider only 
elastically isotropic systems. 



II. THEORY 

We can write the Hamiltonian of the uniform Ising 
model, with allowance for elastic deformations, as fol- 
lows: 
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where S(x) is the n-dimensional order parameter; uo is a 
positive constant; t ~ \T— T c \/T c , where T c is the phase 
transition temperature; u a p is the deformation tensor; 
ai,fl2 are the elastic constants of the crystal; and is 
the quadratic-striction parameter. Passing in equation 
(1) to Fourier transforms and integrating with respect to 
terms depending on non-fluctuating variables, which do 
not interact with the parameter of the S(x) order, and 
introducing for the sake of convenience a new variable 

3 

y( x ) — 2 u aa(x), we obtain the Hamiltonian of the 
system in the following form: 



H = l - J d D q(T + q 2 )S q S_ q 

+ j d D qiS q lS q 2S q 3S-ql- q 2-q3 + 

+a 3 J d D qy ql Sq 2 S- ql -q2 (2) 
a (0) f 1 f la (0) 
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In equation (2) the terms yO, which describe uniform de- 
formations, are separated. As shown in [1], such a sepa- 
ration is necessary, since the nonuniform deformations y q 
are responsible for acoustic phonon exchange and lead to 
long-range effects, which are absent in the case of uniform 
deformations. 

Let us define the effective Hamiltonian of the system, 
which depends only on the strongly fluctuating order pa- 
rameter S, in the following way: 



exp{-H{S]} = B J exp{-H R [S,y}}l[dy q 



(3) 



If an experiment is carried out at constant volume, yo 
is a constant and the integration in equation (3) is only 
done over nonuniform deformations, with uniform defor- 
mations making no contribution to the effective Hamil- 
tonian. At a constant pressure, a term Pfi is added to 
the Hamiltonian, with the volume represented in terms 
of deformation tensor components as 



a=l 



.«W + 0(u 3 )} (4) 



and integration in equaation (3) also performed over uni- 
form deformations. As noted in [6], taking into account 
the quadratic terms in equation (4) can be important in 
the case of high pressures and crystals with large stric- 
tion effects. The neglect of these quadratic terms in [I] 
restricts applicability of the results obtained by Larkin 
and Pikin to only the case of low pressures. As a result, 



we have 

H = \j + q 2 )s q s- g (5) 

d° { qi } S q 1 S q 2 Sq3 S- q 1 _ q 2 - q 3 

+ ^( z o-w ) J d D {qi}S ql S_ ql S q2 S_ q2 , 

z = a?/(4a 3 ), w = a™ 2 / (4a 3 0) ') . 

The effective interaction parameter vq = uq — 12zq, which 
appears in the Hamiltonian owing to the influence of 
striction effects described by the parameter zq, can as- 
sume not only positive, but also negative values. As a 
result, this Hamiltonian describes both first- and second- 
order phase transitions. At v = 0, a tricritical behaviour 
is exhibited by the system. In its turn, the effective in- 
teraction in equation (5), which is determined by the dif- 
ference of the parameters z — w , may lead to a second- 
order phase transition at z — w > and to a first-order 
phase transition at zq — wq < 0. From the given type 
of the effective Hamiltonian follows that there can exist 
a higher-order critical point at which tricritical curves 
intersect if the conditions vq = and zq = wq are ful- 
filled simultaneously [2]. It should be noted that, under 
the tricritical condition zq = wq , the Hamiltonian of the 
model (5) is isomorphic to the Hamiltonian of a uniform 
rigid system. 

In terms of the field-theoretical approach [7], the 
asymptotic critical behaviour and the structure of phase 
diagrams in the fluctuation region are determined by the 
CallanSymanzyk rcnormalization-group equation for the 
vertex parts of irreducible Green function. In order to 
calculate (3- and 7-functions as functions appearing in 
the CallanSymanzyk equation of renormalized interac- 
tion vertices u, a\,a^ or as functions of complex ver- 
tices z — ai/(4a 3 ), w = a^ 2 / (4a^ ) , v = u— 12z, which 
are more convenient for determining the critical and tri- 
critical behaviour of the model, we applied the standard 
method based on Feynmans diagram technique and the 
renormalization procedure [8]. As a result, we obtained 
in terms of the two-loop approximation the following ex- 
pressions for the /^-functions 
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(6) 



n + 2 „ „ 23(n + 2) 2 
-w[ 1 — v - 4nz + 2nw H —— — -v z 



It is known that the perturbation-theory expansions 
are asymptotic and the interaction vertices of order- 
parameter fluctuations in the fluctuation region are suf- 
ficiently large for equations (14) to be applied directly. 
Therefore, in order to extract the necessary physical 
information from the expressions obtained, we use the 
PadeBorel method extended to a three-parameter case. 
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Then the direct and inverse Borel transformations have 
the following form: 

f(v,z,w)= £ a ui2 ,i 3 v h z i2 w i3 

oo 

o 



F(v,z,w)= £ 



{i 1 +i 2 + «3)! 



To obtain an analytical continuation of the Borel trans- 
form of the function, we introduce a series in auxiliary 
variable 9: 



F(v, z, w, 0) 



k=0 



il ,22,23 



fc! 



(7) 

to which we apply the Pade [L/M] approximation at the 
point 8 = 1. This technique was proposed and tested in 
[9] for describing the critical behaviour of a number of 
systems characterized by several interaction vertices of 
order-parameter fluctuations. The property of symme- 
try conservation of the system upon application of the 
Pade approximant in variable 9, revealed in [9], becomes 
significant in describing multivertex models. 

In the two- loop approximation, we used the [2/1] ap- 
proximant for calculating the /3-functions. The nature of 
the critical behaviour is determined by the existence of a 
stable fixed point that satisfies the set of equations 



/3i{v*,z*,w*) = 



(2 = 1,2,3). 



(8) 



The requirement that a fixed point be stable is reduced 
to the condition that the eigenvalues bi of the matrix 



Bi 



d(3j{u\,ul,ul 



(uuUj = v,z,w) (9) 



lie in the right-hand complex half-plane. The fixed point 
with v* = 0, corresponding to the critical behaviour, is 
a saddle point and must be stable in the directions spec- 
ified by the variables z, w and unstable in the direction 
specified by the variable v. The tricritical fixed point is 
stabilized in the direction specified by the variable v as 
a result of taking into account in the effective Hamilto- 
nian of the model terms of sixth order in order-parameter 
fluctuations. The fixed point with z* = w* , which corre- 
sponds to the tricritical behaviour of the second type, is 
also a saddle point and must be stable in the directions 
specified by the variables v, z and unstable in the direc- 
tion specified by the variable w. Its stabilization may be 
due to anharmonic effects. 

The obtained set of the summcd-up functions contains 
a wide variety of fixed points. Table 1 presents fixed 
points for the Ising model (n = 1), XY model (n = 2) 
and the Heisenberg model (n = 3), which are the most 
interesting for describing the critical and tricritical types 
of behaviour, which lie in the physical range of vertices 
with v, z, w > 0. The table also contains the eigenvalues 
of the stability matrix for the corresponding fixed points. 



Analysis of the magnitudes of fixed points and their 
stability suggests the following: the Gaussian fixed points 
10, X0, GO are tricritical and unstable with respect to the 
effect of elastic deformations. The critical behaviour of 
incompressible systems with respect to the deformation 
degrees of freedom is unstable for the Ising model (II) and 
stable for the Heisenberg model (Gl). For theXYmodel 
(XI), the eigenvalue b 2 < 0, but is comparable in order of 
magnitude with the accuracy of calculations; therefore, 
we cannot make any unambiguous conclusion concerning 
the stability of this fixed point. Apparently, the difficul- 
ties in description of the XY model are associated with 
the closeness of the critical dimensionality of the order 
parameter n c to 2. According to the criterion obtained 
in [1], n c < 2, whereas the two-loop approximation gives 
n c = 2.011. For Ising systems, the fixed point is stable at 
constant pressure (12); for Heisenberg systems, the cor- 
responding point is unstable (G2); for the XY model, it 
is impossible to make any unambiguous conclusion since 
the critical dimensionality is also close to 2. The fixed 
points 13, X3 and G3 describe the first type of the tri- 
critical behaviour of compressible systems, which can be 
observed at constant pressure. The fixed points 14, X4 
and G4 are tricritical for systems investigated at constant 
volume. The points 15, X5 and G5 are critical points of 
fourth order, with two critical lines intersecting at these 
points. 

The magnitudes of vertices, obtained in the two-loop 
approximation for the fixed points corresponding to the 
critical and tricritical types of behaviour of the compress- 
ible Ising model, make it possible to calculate the crit- 
ical exponent for the given systems on the basis of the 
expressions (summed by the PadeBorel method) for the 
exponents v and r\: 



HO 



1 /_ n + 2 
12~ 

2(n + 2) , 2 
V = — — — : — V . 



j_ * * ll(n + 2) , 2 

v + nz - nw v 

3888 



243 



(10) 



The values of the other critical exponents can be obtained 
from scaling expressions relating them to the exponents 
v and r\. 

The critical behaviour of compressible Ising systems 
at constant pressure (12) is characterized by renormal- 
ized critical exponents according to Fishers theory of the 
effect of additional thermodynamic variables [10] 



0.632, T] {1) = 0.028, a 



0.103, 



= 0.325, 7 (/) = 1.247. 



For the critical behaviour of the first type (13, X3, G3) 
the Hamiltonian (5) is isomorphic to the Hamiltonian 
of the uniform incompressible model, and, therefore, the 
critical exponents coincide in this case with those for the 
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Table 1. Magnitudes of fixed points and eigenvalues of the stability matrix. 
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V* 


2* 


to* 


bi 


&2 


h 


n=l 


10 











-1 


-1 


-1 


11 


1.064472 








0.6536 


-0.1692 


-0.1692 


12 


1.064472 


0.089187 





0.6536 


0.1702 


0.1710 




1.064472 


0.089187 


0.089187 


0.6536 


0.1702 


-0.1710 


14 





0.5 





-1 


1 


-0.16923 


15 





0.5 


0.5 


-1 


1 


-1 


n=2 


XO 











-1 


-1 


-1 


XI 


0.934982 








0.6673 


-0.0017 


0.1053 


X2 


0.934982 


0.000439 





0.6673 


0.0017 


0.1087 


X3 


0.934982 


0.000439 


0.000439 


0.6673 


0.0017 


-0.1053 


X4 





0.25 





-1 


1 


1 


X5 





0.25 


0.25 


-1 


1 


-1 


n=3 


GO 











-1 


-1 


-1 


Gl 


0.829620 








0.6813 


0.1315 


0.2173 


G2 


0.829620 


0.022909 





0.6813 


-0.1311 


-0.0518 


G3 


0.829620 


0.022909 


0.022909 


0.6813 


-0.1311 


-0.2170 


G4 





1/6 





-1 


1 


1 


G5 





1/6 


1/6 


-1 


1 


-1 



incompressible model: 

v {1) = 0.708, ?7 (7) = 0.028, a (I) = -0.125, 
(3 {I) = 0.364, 7 (/) = 1.397, 
^) = l,^) = 0,a( xy ) = -l, 
/?( xy >=0.5,7 ( * y) = 2, 

/3(G) = o.5, 7 (G) -2. 

The tricritical behaviour of the second type (14, X4, G4) 
corresponds to the critical behaviour of the spherical 
model and is determined by the corresponding exponents: 

^ / ) = l,r ? ( / )=0,a( / ) = -l, 
^) =0 .5,7 (/) = 2, 
^)=l,^) =0l a( xy ) = -l, 
/3< xy )=0.5, 7 (Xy) =2, 



AG) 



(G) 



1. 



l j7 / G ) = o,a 

/3(G) =0.5, 7 ( G > = 2. 



The fixed points of the fourth order (14, X4, G4) are 
characterized by the mean-field values of the critical ex- 
ponents: 

= 0.5, T7 (7) = 0,a (/) = 0.5, 
(3 {I) = 0.25, 7 (/) = 1, 
^ xy )=0.5,r/ xy )=0, a ( xy >=0.5, 
/3^ y )=0.25, 7 (Xy) = l, 
jy( G ) = 0.5,r/ G) = 0,a (G) = 0.5, 

/?( G )=0.25, 7 (G) = 1. 



Thus, the system may exhibit a multicritical behaviour 
under the influence of striction effects. This poses 
the question as to how elastic deformations affect sys- 
tems whose phase diagrams already contain multicritical 
points of bi- or tetracritical nature. Two lines of second- 
order phase transitions and a line of a first-order phase 
transition intersect at the multicritical point in the first 
case, and four lines of second-order phase transitions, in 
the second. In the immediate vicinity of the multicritical 
point, the system manifests a specific critical behaviour 
characterized by the competition of ordering types. In 
this case, one critical parameter is replaced by another 
at the bicritical point, whereas the tetracritical point al- 
lows existence of a mixed phase with coexisting types of 
ordering. Systems of this kind [13] can be described by 
introducing two order parameters transformed in accor- 
dance with different irreducible representations. 

In this case, the model Hamiltonian of the system is as 
follows: 



H = 



+^(*(*)Y + ^(*(*rr + 4 , 

+g iy {x)<S>(x) 2 + g 2 y(x)^(x) 2 + (3y{x) 



(11) 



where $(x) and fy(x) are fluctuating order parameters; 
u w and u 2 o are positive constants; T\ ~ \T — T c i|/T cl , 
r 2 ~ |T — r c2 |/T C 2, where T cl and T c2 are the phase- 
transition temperatures for, respectively, the first and 
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second order parameters; y(x) — u aa{x), where u a p 

a=l 

is the deformation tensor; gi and g 2 are the quadratic- 
striction parameters; (3 is a constant characterizing the 
elastic properties of the crystal; and D is the space di- 
mensionality. As shown above, the striction effects are 
significant only for Ising systems. Therefore, we consider 
here only the case of one-dimensional order parameters. 

Performing transformations analogous to the case of a 
single order parameter, we come to the following effective 
Hamiltonian: 



H =\J dD 1^ + l^^-i 

+ ±Jd D q(T 2 + q 2 )y q y_ q (12) 

-ql-q2-q3) 
f / d D <Zi(* g l* g 2)(* g 3*-,l-,2-,3) 



V02 I aD 
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z\ - w\ 



^ 9i ($ 91 $_ 91 )($ ?2 $_ g2 ) 



+ {z l z 2 -w 1 w 2 ) J d D q t (<P ql <S>- ql )(y q2 y- q2 ) 

voi = uqi - 12z 2 , v 02 = m 2 - 12z|, 
VQ3 = U03 - 12z 1 z 2 , 

9i 92 9i 92 

Z1 = W Z2 = W W1 = W« W2 = Wo 



This Hamiltonian leads to a wide variety of multicritical 
points. As in the case of incompressible systems, both 
tetracritical (v 3 + Yl{ziz 2 — wiw 2 )) 2 < (vi + V2{z\ — 
w\)){v 2 + 12(4 - w 2 

)) and bicritical (v 3 + \2(z\z 2 — 
wiw 2 )) 2 > (vi + \2{z\ - w\)){v 2 + 12(4 - w 2 2 )) types 
of behaviour are possible. In addition, the striction ef- 
fects may lead to multicritical points of higher order. 

In order to calculate functions as functions appearing 
in the CallanSymanzyk equation of renormalized interac- 
tion vertices ui, u 2 , U3, gi, g 2 , <7i°\ g^ or as functions of 
complex vertices zi, z 2 , w\, w 2 , v\, v 2 , v 3 , which are more 
convenient for determining the multicritical behaviour 
of the model, weapplied the standard method based on 
Feynmans diagram technique and on a rcnormalization 
procedure [8]. As a result, we obtained in terms of the 
two-loop approximation the following expressions for the 



functions: 

Pvi = -vi + ^v 2 

Pv2 = ~V 2 + -v\ 



6' 
1 
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6 3 



2 2 1 1 

Pv3 = -Vl + -V 3 + -V1V3 + - 

41 3 23 23 
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23 2 
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8T ? " 


V1V0 — 
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23 o :i 



(13) 
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1 23 , 

0zi = -zi + V1Z1 + 2z\ + 2z 1 z 2 + -v 3 z 2 - gl v i z i 

1 21 . 

^V 2 Z2 



23 2 
W3Z1 - ZT V ' 



Pz2 = ~z 2 + v 2 z 2 + 2z\ + 2z\z 2 + i 

Pwi = -wi + -V1W1 + 4zfwx - 2w\ + Aziz 2 w 2 

-2wiwl + ^v 3 w 2 - - ^HjWi - ^fvjw 2 , 

(3 W 2 = ~w 2 + v 2 w 2 + 4z\w 2 - 2w\ + Aziz 2 wi 
1 23 7 2 

-2w 2 W 2 + -V3W1 ~ -^vlw 2 - 7^ v lw2 ~ 7fj v 3 w l- 



The obtained set of summed-up functions contains a wide 
variety of fixed points, which lie in the physical range of 
vertices with v; t > 0. 

Analysis of the magnitudes of fixed points and of their 
stability suggests the following. The bicritical fixed point 
of incompressible systems (vi = 0.934982, v 2 = 0.934982, 
v 3 = 0.934982, z Y = 0, z 2 = 0, w x = 0, w 2 = 0) is un- 
stable with respect to the effect of uniform deformations 
(61 = 0.090, b 2 = 0.523, b 3 = 0.667, 64 = -0.521, 65 = 
-0.002, b 6 = -0.521, 67 = -0.002). Striction effects sta- 
bilize the tetracritical fixed point for compressible sys- 
tems (vi = 0.934982, v 2 = 0.934982, v 3 = 0.934982, 
Zl = 0, z 2 = 0, wi = 0, w 2 = 0, 61 = 0.090, b 2 = 0.523, 
63 = 0.667, 64 = 2.144, 65 = 0.267, b 6 = 5.223, 
b? = 0.882). 

The question of whether or not other multicritical 
points exist cannot be resolved in terms of the model 
described here since the calculations lead to a degener- 
ate system of equations. The degeneration is lifted when 
terms of higher order in deformation tensor components 
and in fluctuating order parameters are taken into ac- 
count in the Hamiltonian. 

The investigation we performed revealed a significant 
influence of clastic deformations on the critical behaviour 
of compressible systems, which is manifested both in 
a change in the magnitudes of critical exponents for 
the Ising systems and in the appearance of multicriti- 
cal points in the phase diagrams for all the three models. 
The boundary value of the order parameter dimensional- 
ity is close to 2, as in the case of the influence of frozen 
defects. This conclusion is in agreement with the crite- 
rion obtained in [1] on the basis of the mean-field theory. 
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The critical exponents for compressible systems are in 
good agreement with Fishers theory [10] of the influence 
exerted by additional thermodynamic variables. All the 
experimental studies of the multicritical behaviour we are 
aware of have been carried out for compressible systems 
[11]. All of them lead to a conclusion that the critical be- 
haviour is characterized by a Gaussian fixed point with 
mean-field magnitudes of the critical exponents. In a 
number of papers, an opinion has been expressed that the 
system demonstrates a tricritical behaviour [12], while in 



others, it has been stated that this behaviour is bicrit- 
ical. The critical exponents found in this study led us 
to conclude that a tetracritical behaviour determined by 
a fourth-order critical point was observed in the stud- 
ies mentioned above. In systems, described by two fluc- 
tuating order parameters, the striction interaction with 
elastic deformations leads to a change of the bicritical 
behaviour for a tetracritical one. 

The work is supported by Russian Foundation for Ba- 
sic Research N 04-02-16002. 
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